The T=0 neutron-proton pairing correlations in the superdeformed rotational bands 

around 60 Zn 
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The superdeformed bands in 58 Cu, 59 Cu, 60 Zn, and 61 Zn are analyzed within the frameworks of 
the Skyrme-Hartree-Fock as well as Strutinsky-Woods-Saxon total routhian surface methods with 
and without the T=l pairing correlations. It is shown that a consistent description within these 
standard approaches cannot be achieved. A T=0 neutron-proton pairing configuration mixing of 
signature-separated bands in 60 Zn is suggested as a possible solution to the problem. 
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I. INTRODUCTION 

A possibility to find signatures of the T=0 neutron- 
proton (n-p) pairing correlations in N^Z nuclei is re- 
cently a subject of a significant fraction of experimen- 
tal and theoretical studies in nuclear structure physics. 
At low spins, such correlations allow for a consistent de- 
scription of ground states and low-T excitations in even- 
even and odd-odd N=Z nuclei |Q. This type of correla- 
tions may also be, in principle, visible through changes in 
structure of rotational nuclear bands. For example, the 
significance of the so-called delayed alignments in N=Z 
nuclei is at present intensely investigated, both experi- 
mentally d|J|] and theoretically, e.g., see recent Refs. 
and references cited therein. 

In the present study we address another experimental 
fact which may constitute such a signature, namely, an 
anomalous behavior of the second moment of inertia J^" 1 ' 
in the superdeformed (SD) band of 60 Zn, as compared to 
its neighbors. The peak of observed at low spins 
in 60 Zn has been in the original experimental paper |J 
tentatively interpreted as the simultaneous alignment of 
the T=l pairs of gg/ 2 protons and neutrons, although no 
calculation supporting such a hypothesis was presented. 
Together with the discovery of the analogous SD band 
in 61 Zn pQ ], where only a small bump of was ob- 
served, the T—0 paired band crossing was proposed as an 
underlying structure of the 60 Zn band. Indeed, in a sim- 
ple scenario such a crossing would be entirely blocked in 
61 Zn, while for the T—l pairing only the neutron crossing 
would be blocked, leaving half of the peak intact. The 
T=0 paired-band structure was further corroborated by 
the lack of the analogous peak in the SD band in 59 Cu 

On the other hand, the T=l pairing calculations per- 
formed in Ref. E^] indeed resulted in a strong rise of 
with decreasing angular frequency of the 60 Zn SD 
band. However, at lower frequencies solutions could not 
have been obtained, and hence the complete peak of 
was not reproduced. Neither the blocked calculations in 



neighboring odd and odd-odd nuclei were performed to 
support the possibility of reproducing smooth SD bands 
there within the T=\ pairing scenario. It was only ar- 
gued that deformation effects can be important for the 
complete understanding of the physical picture. 

In this study we present the first set of consistent cal- 
culations of the SD bands in 58 Cu, 59 Cu, 60 Zn, and 61 Zn, 
performed within the T=l pairing hypothesis. We show 
that the simple scenario of blocking either the neutron 
or proton T=l pair indeed does not hold, and a more 
complicated picture is obtained. However, a gradual dis- 
appearance of the T=l pairing correlations with increas- 
ing rotational frequency always creates too large values 
of at high frequencies, in disagreement with data. 
In fact at high frequencies the values of , as well as 
the values of relative alignments, are perfectly well de- 
scribed by calculations that altogether neglect the T=l 
pairing correlations. Therefore, it seems that the only 
effect that the no-pairing theory cannot describe is the 
peak of in 60 Zn. Therefore, we attempt to describe 
this structure by a simple T=0 n-p pairing configuration 
mixing of unpaired solutions. 

The superdeformed (SD) bands in the A~60 nuclei 
have already been studied theoretically within various 
approaches Jl3[-pl|. In the present paper we use two 
methods: (i) the cranked Hartree-Fock (HF) method, 
solved by using the HFODD (vl.75r) computer code 
p2fl , with the Skyrme SLy4 p3[ effective interaction 
and no pairing (see Ref. [M for details), and (ii) the 
cranked Strutinsky total routhian surface (TRS) calcula- 
tions based on a deformed Woods-Saxon (WS) potential 
p4j , with the T=l pairing correlations treated within 
the approximate particle number projection by means of 
the Lipkin-Nogami (LN) method (see Refs. &p^| for 
details). Results of these calculations are presented in 
Sees, [n] and III, respectively, while in Sec. [lV| we present 
the T=0 n-p pairing configuration-mixing calculations 
based on the HF results. 
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II. HARTREE-FOCK CALCULATIONS (NO 
PAIRING) 

Before discussing the SD bands in nuclei around 60 Zn, 
we briefly present some generic features of the corre- 
sponding single-particle spectra. The HF neutron single- 
particle orbitals near the SD iV=.Z=30 magic gap, calcu- 
lated in 60 Zn, are shown in Fig. ^. For protons the corre- 
sponding Routhian diagrams are almost identical apart 
from a uniform shift in energy. The single-particle spec- 
tra show large gaps at N(Z)=30 that are stable up to 
the highest frequencies. At the bottom of the SD magic 
gap there appear two strongly deformation-driving in- 
truder orbitals [440]l/2(r=±z), that originate from the 
^0=4 harmonic oscillator (HO) shell, or more specifi- 
cally, from the spherical 1 g 9 / 2 subshell, and therefore are 
denoted as 4 1 = [440]l/2(r=-i) and 4 2 =[440]l/2(r=+i). 
Above the gap, one can see six low-lying orbitals, 
i.e., the next two intruder states 4 3 =[431]3/2(r=— i) 
and 4 4 =[431]3/2(r=-H), as well as four negative-parity 
orbitals which in the present study are denoted as 
/±=[303]7/2(r=±i) and p ± =[310]l/2(r=±i). The /± 
orbitals are in fact the hole states originating from the 
I/7/2 spherical subshell, while the p± orbitals are strong 
mixtures of the 1/ and 2p spherical subshells, i.e., symbol 
p± is assigned only to fix a convenient naming conven- 
tion. 
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FIG. 1. Hartree-Fock neutron single-particle Routhians 
in the SD doubly magic configuration 4 2 4 2 of 60 Zn, calcu- 
lated for the Skyrme interactions SLy4. Lines denoting the 
four (parity, signature) combination are: long-dashed (+,+»), 
solid (+,— i), short-dashed (— ,+i), and dotted (— ,— i). Stan- 
dard Nilsson labels are determined by finding the dominating 
Nilsson components of the HF wave-functions at low (left set) 
and high (right set) rotational frequencies. 

The doubly-magic SD configuration in 60 Zn ||, de- 
noted by 4 2 4 2 , corresponds to occupying all orbitals be- 
low the N—Z—30 gaps, and leaving empty all those that 
are above these gaps. Similarly, following the assign- 
ments of configurations proposed for experimentally ob- 



served bands, we have calculated three other SD bands, 
for the 4H 1 ( 58 Cu @), 4 2 4* ( 59 Cu @), and 4 3 4 2 ( 61 Zn 
JlfH) configurations. The relative alignments (i.e., differ- 
ences of angular momenta at fixed rotational frequencies) 
with respect to the SD band in 58 Cu are shown in Fig. 0. 
Since the experimental SD bands in 59 Cu, 60 Zn, and 61 Zn 
extend to higher rotational frequencies than that in 58 Cu, 
we have artificially extended the latter band by adding 
two gamma rays at 3641 and 4128 keV. This was done 
for the presentation purpose only; alternatively, we could 
have used the 59 Cu band as the reference, however, this 
would have not allowed us to show the relative alignments 
at lower rotational frequencies. Since the exit spins of 
the 58 Cu and 61 Zn bands have been measured only ten- 
tatively, in preparing Fig. we have assumed the values 
of 1=9 and 7=25/2 H, respectively. In calculations, the 
angular momenta I are identified with the average pro- 
jections (Iy). 
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FIG 2. Experimental [|M|u]]l(J and calculated align- 
ments of the SD bands in Cu (diamonds), 60 Zn (squares), 
and 61 Zn (circles), relative to the SD band in 58 Cu (triangles). 
Calculations have been performed within the HF method with 
the SLy4 Skyrme interaction. In order to provide for a suit- 
able reference at high rotational frequencies, two artificial 
gamma rays (open triangles) have been added to the 58 Cu 
experimental data (see text). 

In Fig. U we present a similar comparison between 
the measured and calculated dynamic moments of inertia 
J ( V=dI(uj)/duj. In 58 Cu, 59 Cu, and 61 Zn we obtain very 
good theoretical description of measured relative align- 
ments and second moments. This gives us strong argu- 
ments in favor of the assigned configurations. However, 
unexpectedly, the SD band in the doubly-magic SD nu- 
cleus 60 Zn deviates strongly from the theoretical predic- 
tions. This has been tentatively interpreted as an effect 
of the simultaneous alignment of the g$ / 2 neutrons and 
protons Pi , or as a manifestation of the T=0 n-p corre- 
lations Jul. In the present paper we perform the first 
calculations based on these two assumptions. 
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FIG. 3. Same as in Fig. ^ but for the second moments of 
inertia J^K 



III. STRUTINSKY CALCULATIONS WITH AND 
WITHOUT THE T=l PAIRING CORRELATIONS 

To shed more light especially on the role played by 
the T=l pairing, we performed cranked Strutinsky type 
calculations based on a deformed WS potential. By com- 
paring two sets of calculations, with and without the T=l 
pairing, we aim at tracing the contribution and influence 
of the T=l interaction. The T=l pairing interaction is 
based on a seniority type force and a double stretched 
quadrupole interaction |^] . For the case of odd nucleon 
number and/or excited configurations, each configuration 
is blocked self-consistently pcj . The model has been suc- 
cessful in the description of rotational states in a wide 
range of nuclei. 

To probe the sensitivity of our results to the macro- 
scopic input, we performed two sets of calculations based 
on: (i) the Myers- Swiatecki liquid-drop (MS LD) mass 
formula |n[] and (ii) the folded Yukawa plus exponential 
(FY) mass formula j32|. The MS LD mass formula can 
be considered as rather stiff towards deforming the nu- 
cleus. On the other hand, the FY mass formula, explicitly 
involving the finite range of the nuclear force and the dif- 
fuseness of the nuclear surface, results in a softer surface 
energy and gives larger (3 4 deformations. For very light 
nuclei, the contribution to the surface energy can become 
unphysically large, but for the case of mass ^4=60 region, 
one is still on safe grounds. 

In contrast to the MS LD results, for the FY mass for- 
mula all the four nuclei discussed here have T=l paired 
stable minima at deformations that are comparable to 
those obtained without pairing, but at larger values of 
the hexadecapole deformation parameters. The differ- 
ence in deformations between these nuclei result in dis- 
tinctly different response to the rotating field. Starting 
with 58 Cu, we do not observe any distinct difference be- 
tween the MS LD and FY calculations. Also for the case 
of 60 Zn, no big differences are obtained, although the 
crossing is somewhat sharper here in the FY case. The 



largest difference occurs for the case of 61 Zn. Since the 
FY calculations yield the deformation that is larger than 
for 60 Zn, the neutron g 9 / 2 alignment is becoming more 
smooth, resulting in a rather modest hump in . In 
what follows we concentrate on the results of the FY cal- 
culations. 

In Fig. ^ we present the WS neutron single-particle or- 
bitals near the SD N=Z=30 magic gap in 60 Zn. Even 
though the HF and WS spectra presented in Figs. |l| 
and H have been calculated within so much different ap- 
proaches, they present striking similarities. The equilib- 
rium deformations of the SD shapes, calculated within 
the HF and Strutinsky approaches for the 58 Cu, 59 Cu, 
60 Zn, and 61 Zn nuclei, are presented in Table [| The val- 
ues obtained at hu>=0 and 1 MeV illustrate the degree of 
the rotational polarization occurring along the SD bands. 
Similarly, by comparing the values for the four nuclei one 
can see the effects of the multipole polarizations induced 
by the g g / 2 protons and neutrons, cf. Refs. ||l6| , r2lj| . 
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FIG. 4. Same as in Fig. |l| but for the Woods-Saxon po- 
tential with TRS deformations fa, 7, and fa calculated along 
the T=l paired SD band in 60 Zn. 

The LN method has been shown to be reliable for cal- 
culations of high spin states |^5|,^6|. However, in the 
regime of a very weak pairing, one may encounter nu- 
merical problems in finding a proper solution. Indeed, 
this is the case for the present investigation, where start- 
ing from hui w 1.5 MeV, the static pairing field essentially 
vanishes, and the pairing gaps become of the order of 
100-200 keV. We employ two possible schemes to avoid a 
numerical break down of the paired solution. Either we 
fix the lowest value of the gap parameter to 100-200 keV, 
when no solution is found, or we make a transition to the 
non-pairing calculations. Since the calculations are done 
on a grid in deformation space, the frequency where the 
pairing solution encounters problems differs from point to 
point, giving fluctuations in the total energy. Changes in 
energy of the order of 50 keV are sufficient to cause oscil- 
lations in the calculated moments of inertia. In order to 
address the underlying physics, therefore, we smoothed 
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the moments of inertia in the frequency range where such 
oscillations occur. 

The resulting relative alignments and moments of iner- 
tia are depicted in Figs. || and o, respectively. As can bee 
seen, the strong bump in the experimental moments of 
inertia of 60 Zn is rather well reproduced in the WS+LN 
calculations with T=l pairing. It indeed results from 
the alignment of a pair of <7 9 / 2 protons and neutrons. 
The crossing frequency is somewhat too small in the cal- 
culations, but this can be considered as a detail in this 
context. However, it might also reflect the situation in 
heavier nuclei, where a similar shift has been observed 
H[| and attributed to the lack of the T=0 pairing. Note 
that a similar behavior of J^> has been also obtained 
in the relativistic- mean-field LN calculations of Rcf . jl2| , 
although the increase of J^ 1 ' at low frequencies could not 
have been obtained there. 

TABLE I. Quadrupole {@2 and 7) and hexadecapole (J3a) 
deformation parameters calculated for the SD configurations 
at hu=0 and 1 MeV. For each nucleus the three lines give: 
(a) the HF values obtained from the mass multipole moments 
Q20, Q22, and Qio through the second-order expressions for 



equivalent shapes [|33|34|], (b) the WS potential equilibrium 
deformations obtained by neglecting the pairing correlations, 
and (c) the WS potential equilibrium deformations obtained 
for the LN T=l pairing. In the latter case, the ?iu;=0MeV 
solution in 59 Cu cannot be obtained. 
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58 Cu 


(a) 


0.371 


-1° 


0.051 


0.343 


5° 


0.037 




(b) 


0.392 


0° 


0.038 


0.347 


7° 


0.029 




(c) 


0.374 


-1° 


0.061 


0.357 


6° 


0.024 


59 Cu 


(a) 


0.394 


0° 


0.096 


0.368 


3° 


0.055 




(b) 


0.429 


0° 


0.066 


0.377 


5° 


0.038 




(c) 








0.402 


3° 


0.058 


60 Zn 


(a) 


0.412 


0° 


0.144 


0.391 


2° 


0.089 




(b) 


0.453 


0° 


0.088 


0.418 


3° 


0.058 




(c) 


0.458 


4° 


0.154 


0.426 


2° 


0.089 


61 Zn 


(a) 


0.428 


0° 


0.143 


0.410 


2° 


0.098 




(b) 


0.468 


4° 


0.092 


0.445 


2° 


0.067 




(c) 


0.463 


-1° 


0.123 


0.418 


-1° 


0.085 



correlations due to the Coriolis anti-pairing is of the order 
of 50keV per step in hoj(=0.1) MeV. The sudden drop in 
pairing energy results in a change in the routhian dE^, 
giving rise to this apparent alignment. Thus for all nuclei 
calculated in this study with the T=l pairing, there is 
an excess in the moments of inertia J 1 ^ at high frequen- 
cies, most pronounced in 58 Cu. Such an excess is clearly 
absent in the experimental data. Since 58 Cu is taken as a 
reference for the relative alignment in Fig. ^Jd, the excess 
of alignment obtained in this nucleus for the T=l paired 
calculations perturbs relative alignments shown for other 
nuclei. Another reference nucleus, e.g., 60 Zn may yield a 
better agreement with experiment. 
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FIG. 5. Same as in Fig. ^ but for the Woods-Saxon calcu- 
lations without (upper panel) and with (lower panel) pairing 
correlations. 



The experimental moment of inertia of 58 Cu is totally 
flat, as one would expect since this crossing is blocked. At 
higher frequencies, the calculated J^ 1 moment rises, re- 
sulting in a smaller hump centered at Tiw « 1.5 MeV which 
is absent in the data. In self-consistent calculations, it is 
often difficult to exactly point to the cause of such appar- 
ent alignment as in the case of 58 Cu. The dominant con- 
tribution appears to come from the rather sudden drop in 
pairing energy, where in the region of ?iw=1.4-1.6MeV, 
the pairing gap drops from a value of 0.4 MeV to essen- 
tially zero. At lower frequencies, the change in pairing 



For the case of 61 Zn, the T=l calculations yield min- 
ima in the TRS at large deformation, first after the align- 
ment of the neutron g 9 / 2 orbits. The excess in the mo- 
ments of inertia can again be traced back to the sudden 
drop in the pairing correlations of protons and neutrons. 
Before the alignment of the neutron 39/2, the minimum 
is very shallow at a smaller deformation, where only a 
single 39/2 orbit is occupied. In contrast, calculations 
without pairing yield a minimum that is stable over the 
entire frequency range and has a larger deformation than 
the one in 60 Zn. Finally, the moments of inertia of 59 Cu 
are rather flat, however larger than observed in experi- 
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merit. Again, this is due to the decrease in the pairing 
energy. In 59 Cu, the TRS minimum disappears at low 
frequencies, and therefore the T=l paired band cannot 
be followed to low spins. 
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FIG. 6. Same as in Fig. [] but for the Woods-Saxon calcu- 
lations without (upper panel) and with (lower panel) pairing 
correlations. 

Since we are dealing with nuclei that are located along 
the N—Z line, one may pose the question of the role 
of possible collective T=0 pairing and speculate a little 
about the influence of such correlations. As discussed 
elsewhere |3j|, the collective T=0 pairing field generally 
drives the nucleus to somewhat larger deformation, than 
when only the T=l pair field is present. The sensitiv- 
ity of our results with respect to the macroscopic model 
used, may point to either that the T=l field is too strong 
in our calculations, or by including the T=0 field, the 
results would not be so sensitive to the choice of the ma- 
croscopic model. In addition, since the T=0 pair field 
is more resistant at high angular momenta, one may not 
encounter the unphysical increase in that is present 
in the calculations based on T=l pairing only. 

At low spins, the T=0 field has essentially the same 
properties as the T=l, i.e., resisting the alignment of 
quasi-particles. Assuming that part of the correlations in 
our calculations are indeed due to T=0, would not affect 
much the case of 60 Zn, where we would see a crossing like 
in the calculations with T—l (possibly shifted to some- 



what larger frequencies). However, for the cases of 61 Zn 
and 59 Cu, the blocking effect would be stronger (due to 
the n-p blocking) , and indeed not much of the alignment 
would be observed (as is the case in experiment). 

To really sort out these intriguing problems, unre- 
stricted calculations need to be performed, that simul- 
taneously take into account both T—Q and T=l correla- 
tions. We may however already now conclude, that i) in 
the presence of pairing correlations, one indeed expects a 
hump in the moment of inertia as is observed for the case 
of 60 Zn and ii) the simple blocking picture does not hold 
here, where strong polarizing effects are present, yielding 
different deformation for the nuclei discussed here and 
as a result, different pattern of the alignment. Before 
such complete solutions become available, and the ex- 
pectations expressed above can be corroborated, in the 
next section we investigate a very simple non-collective 
T=0 n-p pairing scenario by considering the configura- 
tion mixing of unpaired HF solutions. 



IV. HARTREE-FOCK CALCULATIONS WITH 
THE T=Q N-P PAIRING CONFIGURATION 
MIXING 

Apart from the 4 2 4 2 configuration discussed above, in 
60 Zn we also calculated 6 other configurations, namely 
those that correspond to exciting the 4 2 proton and neu- 
tron simultaneously to the negative-parity orbitals f± 
and p±. In principle, there are 16 such excitations pos- 
sible, however, the lowest ones are obtained by putting 
the neutron and the proton in the same orbitals. This 
gives 4 configurations denoted by 4 1 / + 4 1 / + , 4 1 /_4 1 /_, 
4 1 p + 4 1 p + , and 4 1 p_4 1 p_. In addition, we also study 2 
other configurations obtained by putting the neutron and 
the proton into the [303] 7/2 orbital with different signa- 
tures, i.e., those denoted by 4 1 / + 4 1 /_ and 4 1 /_4 1 / + . 

In Fig. |^(b) the energies of the seven configurations 
selected above are shown with respect to a common 
rigid-rotor reference energy of 0.025x7(7-1-1) MeV. Sim- 
ilarly, Figs. 0(a) and (c) show the analogous configura- 
tions in 58 Cu and 62 Ga. Because orbitals f± and p± 
are very close in energy (cf. Fig. |l|), they strongly inter- 
act and mix, which very often precludes the convergence 
of the HF procedure, see discussion in Ref. |56|. Apart 
from that, the bands of Fig. ^ are shown up to the so- 
called termination points, i.e., up to the point where the 
angular-momentum contents of the involved orbitals does 
not allow for a further angular momentum build up, see 
Ref. |^?| , without a significant rearrangement of the nu- 
cleons. 

By considering the available projections of the total 
angular momentum I y for oblate shapes with the y axis 
as the symmetry axis, one can easily determine the values 
of the termination-point angular momenta 7 t , see Table 
|j. The bands obtained in the HF calculations do not 
always terminate at the oblate axis and can usually be 
continued beyond 7 t . However, at angular momenta 7 t 
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there always occur significant changes in the structure of 
bands. Below we discuss and present results only up to 
the termination points I t . 
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FIG. 7. Energies of selected configurations in 58 Cu (a), 
60 Zn (b), and B2 Ga (c), calculated within the HF method 
with the Skyrme Shy4 interaction, and plotted with respect 
to a rigid-rotor reference energy. The configurations shown 
in the legend correspond to the numbers of occupied iVo=4 
intruder orbitals, n and p, that are indicated in each panel. 

A conspicuous feature of the HF energies presented in 
Fig. [7] is the significant energy separation between the 
n-p paired configurations 4™/+4 p /+ and 4™/_4 ?) /_ on 
one side, and the broken-pair configurations 4™/+4 p /_ 
and 4 n /_4 p / + on the other side. The former and lat- 
ter configurations have opposite total signatures, i.e., in 
the even-even nucleus 60 Zn, configurations 4 1 / + 4 1 / + and 
4 1 i /L4 1 /_ (4 1 /+4 1 /- and4 1 /-4 1 /+) correspond to r=+l 
(r=— 1), while in the odd-odd nuclei 58 Cu and 62 Ga the 
analogous configurations correspond to r=— 1 (r=+l). 
Such a signature-separation effect has been for the first 
time discussed for the SD bands in 32 S |3S[| . Here it is 
obtained in the heavier SD region of the A~60 nuclei, as 
a mutatis mutandis identical effect occurring for all the 



orbitals promoted to the next HO shell. 

TABLE II. Values of the termination-point angular mo- 
menta It (in h) for the seven selected configurations in 58 Cu, 
60 Zn, and 62 Ga. For convenience, the second column gives 
the configurations shown in the convention of Refs. J3^|,[L7|, 
that, however, does not allow for distinguishing between the 
signatures of the occupied orbitals. 



Configuration 


58 Cu 
n=p=0 


60 Zn 
n—p=l 


62 Ga 
n=p—2 




[2(p+l),2(n+l)] 


29 


36 


41 


4"/+4 p /+ 


[lp, In] 


15 


24 


31 


4 n /„4 p /_ 


[lp, In] 


13 


22 


29 


4 n p+4?p+ 


[2p, 2n] 


23 


32 


39 


4>_4 p p_ 


[2p, 2n] 


21 


30 


37 


4 n / + 4 J 7_ 


[lp, In] 


14 


23 


30 


4'7-4 P /+ 


[lp, In] 


14 


23 


30 



In Ref. |g8j the signature-separation effect was inter- 
preted as a result of the strong n-p attraction transmitted 
through the time-odd mean fields. Such an attraction is 
typical for any realistic effective interaction, and it has 
its origin in the spin-spin components of the interaction. 
(The signature separation vanishes when in the Skyrme 
energy functional |]39|] the coupling constants correspond- 
ing to terms s-s and s-As are set equal to zero.) When 
averaged within the mean-field approximation, the spin- 
spin c omp onents lead naturally to the time-odd mean 
fields p9| . Within the phenomenological mean fields, 
like those given by the Woods-Saxon or Nilsson poten- 
tials , the time-odd mean fields vanish, and therefore 
all the four configurations 4 ra /±4 p /± are nearly degener- 
ate, i.e., the signature-separation effect occurs only for 
self-consistent mean fields generated from the spin-spin 
interactions. 

One should note that the four configurations 4™/±4 p /± 
have purely independent-particle character (Slater- 
determinant wave functions), i.e., no collective pair cor- 
relations are built into the wave functions. Neverthe- 
less, configurations 4™/ + 4 p / + and 4 n /_4 p /_ contain one 
more T=0 n-p pair as compared to the 4"/ + 4 p /_ and 
4 n /„4 p / + configurations, and therefore are sensitive to 
the n-p pairing component of the effective interaction 
that is attractive. As a result, the paired configurations 
4 n / + 4 p / + and 4 n /_4 p /_ cross the magic configurations 

4 n+l 4 p+l at J=llj lg) and 27 fr i n 58 CU) 60 Zn) and 62 G&j 

respectively. 

The n-p pairing correlations should be, in principle, 
studied by using methods beyond the mean-field approx- 
imation, i.e., by taking into account the configuration- 
mixing effects for configurations that differ by the n- 
p pair occupations. The generator-coordinate method 
(GCM) is the approach of choice for including such 
effects. It allows for a consistent improvement of wave 
functions, while staying in the framework of the vari- 



G 



ational approach. Therefore, the same interaction can 
be/should be used in the HF method and in the mixing 
of the HF configurations via the GCM method. 

At present, the GCM approach in the rotating frame 
has not yet been implemented, and in the present study 
we discuss the same physics problem by introducing a 
model T=0 n-p pair-interaction Hamiltonian in the form 
of 



^n-p -^0 H~ ^n-p ^ ^rar-^rar ^ 

rar ct/3r 



0T1 



(1) 



where the particle-number (N Tar ) and T=0 n-p pair- 
creation (Pl r ) operators read 



Nrar 



t car ■ t "Otr 1 



(2a) 
(2b) 



t denotes neutrons (y) or protons (7r), and a and /3 de- 
note the Nilsson labels without the signature quantum 
number r that is shown explicitly. 

Hamiltonian (Q) is meant to replace the usual effective- 
interaction (Skyrme) Hamiltonian when studying the n- 
p correlation aspects of the nuclear wave functions, and 
not to be added on top of it. Therefore, the effective 
single-particle energies e rQr and the coupling constants 
G"^ p =G^ a p have to be angular-momentum and config- 
uration dependent, and Hamiltonian ([!]) should be un- 
derstood as a phenomenological interaction operator be- 
tween configurations that differ by the n-p pair occupa- 
tions. 

The diagonal pairing term can be transformed as 



pt p 

ar ar 



N l , ar N 7Tar , 



(3) 



i.e., it gives a non-zero contribution only if both a neutron 
and a proton occupy the given {ar} orbital. Therefore, 
the diagonal matrix elements of the n-p pairing Hamilto- 
nian (pi) in any given configuration, 



E^conf.) = (conf.|i/ n _ p |conf.] 



(4) 



can be immediately calculated for each Slater determi- 
nant. In particular, since in all the four configurations 
4 n /±4 p /± the effective single-particle energies are iden- 
tical (cf. the Routhian diagram in Fig. [I]) , the differences 
of the total energies in the signature-separated configu- 
rations read 



£(4™/ + 4*>/_) - E(4"/-4 P /+) = 0, 
£(4"/+4P/+) - B(4"/_4f/_) = 0, 
E(4"/ + 4f/_) - £(4"/+4P/+) = G n _p[303]7/2, 



(5a) 
(5b) 
(5c) 

where G n - P [303]7/2 stands for the diagonal matrix ele- 
ment G« Q p for a= [303] 7/2. 

By subtracting the total HF energies of configurations 
in Eq. (|kj), see Fig. [?], one thus obtains an estimate of 
the n-p pairing diagonal matrix element G^°L. Such rel- 
ative energies ( |5c| ) in 58 Cu, 60 Zn, and 62 Ga are plotted 



in Fig. H One can see that the effective matrix elements 
depend strongly on the angular momentum, and decrease 
from G^ p (I=0)^1.6 ( 58 Cu) or 1.9 MeV ( 60 Zn and 62 Ga), 
reaching zero at the termination angular momentum I t . 
This dependence can be very well parameterized by a 
simple cubic expression, 



G™(I) = G£»(J=0) x 



shown by dashed lines in Fig. 



2.0 




(6) 



> 

^ 1.5 

>. 

O) 

k. 

a 

c 1.0 

LU 

> 
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FIG. 8. Energies of configurations 4™/+4 p /_ (closed cir- 
cles) and 4 n /_4 p /+ (open circles) in 58 Cu, 60 Zn, and 
62 Ga, relative to the corresponding 4"'/+4 p /+ configura- 
tions. Dashed lines show the simple cubic approximations 
of Eq. The relative energies can be identified with the 

angular-momentum-dependent T—0 n-p pairing matrix ele- 
ments G n -p in the [303] 7/2 orbital (see text). 

A large standard signature splitting of the other single- 
particle orbitals, which have lower values of the K quan- 
tum numbers, does not allow us to determine the other 
diagonal matrix elements G^p directly from the HF re- 
sults, as in Eq. (fcj). Of course, such a determination 
of the non-diagonal matrix elements is not possible ei- 
ther. However, we may use the /-dependence of Eq. (|^) 
to postulate a simple separable approximation for the n-p 
pairing interaction matrix G%?L in the form 



Gf p (I) = y/G^(I)G^(I). 



(J) 



Such a postulate is motivated by the fact that the pair- 
ing matrix elements of short-range interactions are given 
primarily by overlaps between the space wave functions, 
or more precisely, by the integrals of products of squares 
of the wave functions. Then, Eq. (|7]) stems from ap- 
proximating the integral of products by the product of 
integrals. 

Within the separable approximation (0), the T=0 n- 
p pairing interaction V n . p in Hamiltonian ([!]) takes the 
simple form of 



7 



V,, 



-G(J)(p^P+i + (8) 
where the /-dependent collective n-p pair operators read, 
P} = J2*M)Pl, (9) 



for 



x 2 a (I) = G™(I)/G(I), 
G(I) =Y,GT P (I)- 



(10a) 
(10b) 



Even then, however, the problem is defined by one pa- 
rameter per orbital, G"_ p (7=0), i.e., it cannot be defined 
without explicit microscopic configuration-mixing calcu- 
lations. 



-502 



-504 



-506 




10 20 30 

Angular momentum 1(h) 

FIG. 9. Same as in Fig. (?](b) but for configurations 
4 2 4 2 and 4 1 /+4 1 /\j_ interacting through the T=0 n-p pair- 
ing interaction || with G n - P [303]7/2(/=0) = 1.9 MeV and 
G n -p[440]l/2(7=0)=0.65MeV. Crosses show the experimen- 
tal data in the absolute energy scale. The symbols used in the 
figure indicate which are the pure configurations of Fig. |^(b) 
that dominate in the given mixed configurations. 
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FIG. 10. Experimental (full squares) and calculated (open 
squares) dynamic moments of inertia in the SD band of 
60 Zn. Calculations correspond to the lowest band shown in 
Fig- | 



Before these become available, in the present study 
we perform the simplest two-level mixing calculation, in 
which the two configurations that cross in 60 Zn, 4 2 4 2 and 
4 1 / + 4 1 / + , see Fig. 0(b), are allowed to interact through 
the T—0 n-p pairing interaction (B). With the diagonal 
matrix elements of Hamiltonian (1) taken from the HF 
calculations, and the interaction matrix element defined 
by the value of G n _ p [303]7/2(7=0) = 1.9MeV, also taken 
from the HF calculations, we are left with one free pa- 
rameter, i.e., with the value of G n _ p [440]l/2(7=0). 

By fixing this parameter at G n - P [440] 1/2(7=0) = 
0.65 MeV, we obtain at the crossing point of 1=18 h the 
effective interaction strength of 0.79 MeV. With the I- 
dependent matrix elements given by Eqs. (||) and (0), 
we obtain the energies and dynamic moments of inertia 
shown in Figs. |t] and |l0|, respectively. It is clear that the 
mixing and interaction of the 4 2 4 2 and 4 1 / + 4 1 / + config- 
urations correctly reproduces the magnitude of the bump 
in the J {2) of 60 Zn. 

The position of the crossing point is obtained at fre- 
quency or spin that are too large by 0.2 MeV or Ah, 
respectively, as compared to experiment. As seen in 
Fig. 0(b), this position is dictated by the diagonal ma- 
trix element G n - P [303]7/2 that shifts down configuration 
4 1 / + 4 1 / + with respect to the broken-pair degenerate con- 
figurations 4 1 / + 4 1 /_ and 4 1 /_4 1 / + . As discussed above, 
such a shift is a direct consequence of the time-odd mean 
fields resulting from the Skyrme energy density. In the 
present work we have used the time-odd terms as directly 
given by the SLy4 Skyrme functional, see Ref. [|3Sf| , i.e., 
those that result from fitting the time-even, and not time- 
odd properties of nuclei. It is clear that a modification 
of these time-odd terms, that is permitted in the local 
density approximation, may move the crossing frequency 
from its current position in Fig. |l^. In fact, it is ob- 
vious that by decreasing this intensity one may easily 
decrease the crossing frequency. We do not attempt such 
a fit here, because the problem of finding good physical 
values of the time-odd coupling constants is much more 
general, and it would not make too much sense to make 
such an adjustment based solely on the specific effect 
discussed in the present study. We only note in passing 
that an analogous readjustment of the isovector time-odd 
coupling constants Q has led to values that are quite 
different from those resulting directly from the Skyrme 
functional. 



V. CONCLUSIONS 

In the present study we presented a model description 
of the four SD bands in nuclei around 60 Zn in order to 
evaluate the role of the T=l and T=0 pairing correla- 
tions at high spin in N~Z nuclei. On the one hand, we 
have shown that calculations with no pairing, whether 
within the Strutinsky- Woods-Saxon or Skyrme-Hartree- 
Fock approaches, provide an excellent description of all 



8 



bands except the one in the doubly-magic N=Z nucleus 
60 Zn. On the other hand, the standard Lipkin-Nogami 
treatment of the T=l pairing gives a fair description of 
the bump in the second moment of inertia in 60 Zn, 
which results from the simultaneous alignment of the go, 12 
pairs of neutrons and protons, but fails in describing low 
values of in all the four nuclei at high frequencies. 
The latter effect results from a gradual disappearance of 
the T—l pairing correlations with increasing spin, and 
cannot be avoided if these same pairing correlations have 
to be responsible for the positive result in 60 Zn. We have 
also shown that the deformation changes caused by polar- 
ization effects of high-j orbitals are strong, and strongly 
modify the simple blocking picture when going from even 
to odd isotopes. Nevertheless, even with these polariza- 
tion effects taken into account in a self-consistent way, 
the overall description of the discussed set of bands is 
not satisfactory. 

In looking for an alternative physical scenario we have 
shown that another kind of crossing results from the 
signature-separation effect that shifts down the HF con- 
figurations with neutron-proton pairs present, with re- 
spect to those were such pairs are broken. Then the 
n-p paired configurations cross and interact, giving a cor- 
rect qualitative reproduction of in 60 Zn. This sce- 
nario has the advantage that the effect of interaction en- 
tirely disappears whenever an odd neutron or an odd pro- 
ton blocks the n-p pairing interaction. Therefore, non- 
interacting configurations are obtained in odd neighbors 
of 60 Zn, yielding a perfect agreement with the data. This 
result was, however, obtained within a very simple two- 
band interaction and two-band mixing scheme, while the 
collective features of the T=Q n-p pairing were not stud- 
ied. Whether or not such collective aspects of a simulta- 
neous presence of the T—l and T—0 pairing correlations 
are important will be the subject of future investigations. 
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